Anderson localization of entangled photons in an integrated quantum walk 
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Waves fail to propagate in random media. First 
predicted for quantum particles in the presence of 
a disordered potential, Anderson localization has 
been observed also in classical acoustics, electro- 
magnetism and optics. Here, for the first time, we 
report the observation of Anderson localization 
of pairs of entangled photons in a two-particle 
discrete quantum walk affected by position de- 
pendent disorder. A quantum walk on a disor- 
dered lattice is realized by an integrated array of 
interferometers fabricated in glass by femtosec- 
ond laser writing. A novel technique is used to 
introduce a controlled phase shift into each unit 
mesh of the network. Polarization entanglement 
is exploited to simulate the different symmetries 
of the two-walker system. We are thus able to 
experimentally investigate the genuine effect of 
(bosonic and fermionic) statistics in the absence 
of interaction between the particles. We will show 
how different types of randomness and the sym- 
metry of the wave-function affect the localization 
of the entangled walkers. 

In 1958 P.W. Anderson [1 predicted that the wave- 
function of a quantum particle can be localized in the 
presence of a static disordered potential. As a conse- 
quence of this mechanism it is expected that particle and 
energy transport through a disordered medium should be 
strongly suppressed and that an initially localized wave 
packet should not spread out with time. After more than 
fifty years from its discovery Anderson localization is still 
widely studied and it has pervaded many different ar- 
eas of physics ranging from condensed matter and cold 
atomic physics to wave dynamics and quantum chaos [2] . 
This phenomenon emerges quite generically in the behav- 
ior of waves in complex media, and it has been experi- 
mentally observed in a variety of different systems: Bose- 
Einsten condensates 01 14] , light in semiconductor pow- 
ders [5l and photonic lattices [6HH] , single photons in bulk 
optics [9 and in fiber loops [10], microwaves in strongly 
scattering samples [11], besides ultrasound waves in a 
three-dimensional elastic system [12j. 
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FIG. 1: Concept scheme of cascaded beam splitters to im- 
plement photonic quantum walks. Disorder is introduced by 
phase shifters placed at each beam splitter's output port, be- 
fore entering the next one. 



Anderson localization is a single-particle process which 
arises from the destructive interference among different 
scattering paths. Nevertheless, even in the absence of a 
direct interaction between particles, pure quantum cor- 
relations [13^ are expected to influence in a non-trivial 
way the underlying localization dynamics p!4HT7] . By 
taking advantage of the perfect phase stability provided 
by miniaturized integrated waveguide circuits [18 , we ex- 
perimentally simulate a quantum walk of a two-photon 
polarization-entangled state in a disordered medium. We 
are thus able, through a mapping derived in Ref. [14 , to 
test the localization of a pair of non interacting particles 
obeying bosonic/fermionic statistics [19 . 

A quantum walk (QW) [20 is an extension of the clas- 
sical random walk, where the walker goes back and forth 
along a line and the direction at each step depends on the 
result of a fair coin flip. At the quantum level, interfer- 
ence and superposition phenomena lead to a non-classical 
behavior of the walker giving rise to new interesting ef- 
fects, which can be harnessed to exponentially speed up 
search algorithms [21 and to realize universal quantum 
computation [22 . Besides, QWs have also been proposed 
to analyze energy transport in biological systems 231 El] • 
Different experimental implementations of single-particle 
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FIG. 2: (a) Scheme of the network of directional couplers implementing a 8-step quantum walk with static disorder. Different 
colors stand for different phase shifts, (b) Controlled deformation of either of the two S-bent waveguides at the output of each 
directional coupler extends the optical path and is equivalent to the application of a phase shifter, (c) The deformation is given 
by a non-linear coordinate transformation, which is function of a deformation coefficient d (see Supplementary Information). 
The graph shows the undeformed S-bend (solid line), together with a deformed one (dashed), (d) Schematic of the Mach-Zehnder 
structure, representing the unit cell of the directional couplers network, fabricated for calibrating the phase shift induced by the 
deformation, (e) Phase shift induced by the deformation: theoretical curve calculated from the nominal geometric deformation 
(solid line), and experimental measurements (diamonds). 



quantum walks were performed with trapped atoms (2^ , 
ions [26l [27] , energy levels in NMR schemes [28] , photons 
in waveguide structures [29j , in bulk optics [9l [30] , and 
in a fiber loop configuration [lOl [HI [32] . Very recently 
quantum walks of two identical photons have been per- 
formed only in ordered structures [I8l [33l [34] . 
A physical realization of a discrete QW can be provided 
by photons passing through a cascade of balanced beam 
splitters (BSs) arranged in a network of Mach-Zehnder 
(MZ) interferometers as shown in Figure [l] Here each 
BS implements simultaneously the quantum coin oper- 
ation, i.e. the choice of the direction the particle will 
move in, and the step operator, which shifts the walker 
in the direction fixed by the quantum coin state (the time 
evolution being simulated stroboscopically) [18 . Ac- 
cordingly every output of a BS of the network corre- 
sponds to a given point in the space-time of the QW, 
the horizontal rows of the setup representing different 
time steps. In this scenario, disorder can be added in the 
QW evolution by simply introducing (randomly selected) 
phase shifts between the MZ interferometers paths (see 
Fig.jl]). In particular the time-independent, static^ disor- 
der needed to enforce Anderson localization on the pho- 
tonic walker, is obtained by fixing the same phases for 
all the MZs which correspond to the same lattice site. 
Making sure that response of the device is polarization 
independent, the localization of a non interacting entan- 
gled pair can now be studied by injecting in two different 
ports of the device -namely A and B in Fig [2^- a two- 
photon entangled-polarized state generated via sponta- 
neous parametric down- conversion. The above approach 



would be extremely hard to implement with bulk optics 
mainly because of size and of very challenging stability 
issues. However in the last few years integrated quan- 
tum photonics proved to be a highly promising exper- 
imental platform for quantum information science [36]. 
Recently, integrated waveguide circuits have been em- 
ployed for quantum applications, in order to realize two- 
qubit gates [35j [37H39] , quantum algorithms [40 , quan- 
tum walk on a chip [18l [33l [41] and enhanced quantum 
sensitivity in phase-controlled interferometers J42H44] . In 
order to observe Anderson localization for polarization 
entangled photons an important step forward is required 
in the available experimental platform. Our setup, for 
the first time, integrate all the necessary ingredients to 
this aim: polarization independent elements, interfero- 
metric structures, together with a proven capability to 
implement suitable phase shift in different points of the 
QW circuit. Furthermore, in order to get a convincing 
evidence of localization it has been necessary to observe 
the experimental simulation for quantum walk of differ- 
ent steps. We were able to realize up to eight-step QW 
circuits affected by a controlled disorder, thus integrating 
tens of BSs on the same chip. 

In our experiment the setup of Figure [l] has been re- 
alized by using integrated waveguide circuits, as in Fig- 
ure |2^, where BS elements are replaced by directional 
couplers. The discrete Tn-axis indicates the different sites 
of the QW, while the discrete n-axis identifies the differ- 
ent time steps. The integrated waveguide circuits have 
been fabricated by femtosecond laser writing technol- 
ogy ^45, ,464 . This technology exploits nonlinear absorp- 



tion of femtosecond pulses, focused below the surface of 
a transparent dielectric substrate, to obtain a perma- 
nent and localized refractive index increase. Transla- 
tion of the sample under the laser beam along the de- 
sired path enables the fabrication of optical waveguide 
circuits with arbitrary three-dimensional geometries. In 
order to obtain a totally polarization independent be- 
havior, the 3-dimensional geometry detailed in Ref. [18] 
has been adopted (these devices are known to allow the 
propagation of polarization entangled states [35 ). The 
phase shifters are implemented by deforming one of the 
S-bent waveguides at the output of each directional cou- 
pler (green or red segments in Figure [2J3) , in order to 
stretch the optical path. The phase shift [— tt, tt] in each 
MZ cell is implemented by lengthening the optical path in 
the green segment to introduce a [0,7r] phase shift, while 
the complementary range [— tt, 0] is achieved by lengthen- 
ing the red segment. In this way, smaller deformations, 
always of the same kind (lengthening of the path), are 
capable to provide the full range of phase shifts. 

Figure ^ shows both an undeformed and a deformed 
S-bend. 
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Steps 


Input A 


Input B 
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0.991 ±0.002 


0.992 ± 0.002 


6 


0.994 ±0.003 


0.987 ±0.003 


8 


0.951 ±0.004 


0.946 ± 0.005 


STATIC 


Steps 


Input A 


Input B 


4 


0.980 ±0.003 


0.976 ± 0.002 


6 


0.985 ± 0.002 


0.976 ± 0.003 


8 


0.938 ± 0.004 


0.957 ±0.004 



TABLE I: Similarities for single particle QW distributions. 
The values are calculated as mean average on distributions of 
single-photon in different polarization states. 



To test our technique and calibrate the achieved phase 
shift as a function of the imposed deformation d (see Sup- 
plementary Information for a detailed definition) several 
MZ interferometers were fabricated with the design of 
Figure 2] d) , reproducing exactly the unit cell of the QW 
network. Each interferometer has one S-bend (the one 
colored in the figure) deformed with a different value of 
d. Laser light at A =806 nm wavelength was injected 
in the interferometers and the induced phase shift was 
then retrieved from the measured light distribution at 
the output. Figure [2] e) reports the experimentally mea- 
sured phase shifts as a function of the deformation pa- 
rameter d. The experimental points are in good agree- 
ment with the phase shift predicted by evaluating nu- 
merically the geometric lengthening A/ of the deformed 
S-bend 4>theo = x^^- 

We implemented a lattice with static disorder by im- 
posing the same phase shift to the MZ cells correspond- 



ing to a fixed site of the QW line as in Fig. |2fe 
(prni^^)' QW circuits composed by 4, 6 and 8 step af- 
fected by static disorder were realized in a way that the 
4-step phase pattern was embedded within the 6-step 
phase pattern and, in turn, this was embedded within 
the 8-step one. 

Another set of 4, 6 and 8 step ordered, i.e. with per- 
fectly symmetric MZ cells, QW circuits was realized and 
compared with the corresponding disordered one. 

First of all we measured the single particle distri- 
butions (see Supplementary Information) in order to 
demonstrate the polarization insensitivity of the inte- 
grated QWs. We repeated this measurement by in- 
jecting single photon states with different polarizations. 
The measured distributions exhibit very similar behav- 
iors. We compared the obtained results with the ex- 
pected ones by calculating the similarity defined as 5 = 



(Ei,i y Aj^-,^) V Ei,i Aj Ei,i D[., which is a general- 
ization of the classical fidelity between two distributions 
D and D' . The obtained values are reported in Table 
|l] for the ordered QW circuits and for QWs with static 
disorder. These high values and low deviations highlight 
once more the fabrication control and polarization insen- 
sitivity of our integrated devices. 

As a second step we carried out the investigation of 
two-particle QWs. The investigation of Anderson lo- 
calization for bosonic and fermionic particles was real- 
ized by adopting the complete experimental setup de- 
scribed in details in the Supplementary Information. As 
mentioned above, we exploited polarization entanglement 
to simulate bosonic and fermionic statistics. To this 
aim polarization-entangled photon pairs, generated via 
spontaneous parametric down-conversion, were simulta- 
neously injected into arms A and B of the 4, 6 and 8 steps 
QW circuits to observe the progressive quench of photon 
propagation in disordered QWs (Fig. |2k). By setting 
the phase (/) of the state -^[\H)a\V)b ± e^^lV)^!^^)^], 
bosonic ((/) = 0) or fermionic ((/) = tt) QWs were observed. 
In Figure|3]we show how entangled photons localize in the 
presence of a random static potential by plotting the joint 
probability Pj^k of detecting one particle in the output 
port j and the other in the output port k (the probabil- 
ity being determined by collecting events independently 
from the photon polarization) . The different panels com- 
pare the ordered and disordered cases in the case of sym- 
metric and antisymmetric wave-functions. We report 
also the case of single photons which are reconstructed by 
tracing out the position of one of the particles of the en- 
tangled pair. While in the case of an ordered system the 
walkers spread on with the increasing number of steps, 
Anderson localization implies that the wave-packets will 
remain localized around the central sites irrespectively of 
the number of steps. This is indeed what we observed. 
The difference between the ordered and disordered case is 
most evident for the 8-step QW, compare the panels (g) 
and (h) with (r) and (s) in Figure [3] The agreement of 
the experimental data with the theoretical predictions. 
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FIG. 3: Experimental results of single- and two-photon distributions for bosons and fermions in an ordered QW (a-i) and in a 
QW in presence of static disorder (1-t). Single-particle density distributions have been computed by tracing out the position of 
one of the particles (summing over the columns of the two-photon probability distribution matrices). 



ORDERED 


STATIC 


Steps 


Bosons 


Fermions 


Single 


Steps 


Bosons 


Fermions 


Single 


4 


0.946 ± 0.003 


0.914 ±0.003 


0.996 ±0.001 


4 


0.918 ±0.003 


0.902 ±0.003 


0.993 ±0.001 


6 


0.940 =b 0.003 


0.851 ±0.003 


0.997 ±0.001 


6 


0.890 ±0.006 


0.903 ± 0.004 


0.985 ± 0.002 


8 


0.768 ± 0.006 


0.780 ± 0.007 


0.934 ± 0.004 


8 


0.803 ± 0.004 


0.785 ± 0.004 


0.947 ±0.002 



TABLE II: Similarities between the experimental distributions of Fig. [3] and the expected ones. Uncertainties arise from the 
Poisson distribution of counting statistics. 



again quantified by the similarities, is reported in Ta- 
ble [nj In the ordered case S is slightly worse for the 
8-step QW. This discrepancy, due to some unavoidable 
uncertainty in BSs realization is milder in the disordered 
case. Here, as expected, additional phase-shift to the 
"intentionally-chosen" random one will have less effect 
due to localization. 

The entangled pairs localize in a manner which de- 
pends on their statistics. A more quantitative estimate 
of the difference in the localization properties of entan- 
gled photons may be obtained by looking at the variance 
of the two-photon mean position xm = {j ± k)/2 associ- 
ated with the probability distributions of Figure [3] as a 
function of the number of steps. This is shown in Fig- 
ure [4k, b. While for the ordered case the variance grows 



quadratically with the number of steps n, it is weakly de- 
pendent on n in the disordered case, indicating that the 
system tends towards localization. The numerical simu- 
lation, performed by considering discrete time QWs with 
static disorder shows that localization starts even with 
a relatively small number of steps, giving rise to a fully 
localized state (corresponding to a variance not varying 
with n) typically for n ~ 100. However, such number 
of steps is currently out-of-reach for any technological 
platform. The clear difference observed in our experi- 
ment between ordered and disordered QWs is a strong 
evidence of the onset of Anderson localization for a small 
number of steps. 

In addition we observe a very interesting feature emerg- 
ing from the data, both in Fig. [3] and in Fig. |4| further 



(a) 

— 10 

s 

I 



(c) 



• bS 


■ f:: 


Theo 


Theo 
"^Ord 


o b:: 


□c 



50 



25 






■-Br 

Theo 
I" Stat 


02 

(0 

1^ 


• 

■ ^ - ' 



4 6 8 

Number of steps 



4 6 

Number of steps 



4 6 8 

Number of steps 




Bosons Ordered ,, 



H 



(f) 



Fermions Ordered 



llHililii 



0123456789 10 11 12 13 14 15 01234567 



12 13 14 15 



(g) 



I Bosons Static 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

R R 



FIG. 4: (a-b) Variance of the two photon mean position 
xm = (j + k)/2 and (c,d) of the relative distance R = K ~ ^| 
shown in function of the number of steps. Experimental re- 
sults are reported for bosons (circles) and fermions (squares). 
Empty markers refer to ordered QWs, whereas filled ones cor- 
respond to QWs with static disorder. In panels a,b dot-dashed 
and solid lines represent theoretical behaviors for bosons and 
fermions in the ordered case while in panels c,d dotted and 
dashed lines correspond to theory for bosons and fermions 
in presence of static disorder, respectively. Error bars are 
smaller than the spot size, (e-h) Probability distributions of 
the relative distance R for bosonic (left) and fermionic (right) 
two-photon 8-step QWs in the case of a ordered QW (e-f ) and 
in presence of position-dependent disorder (g-h) . 



evidencing the different behaviour between bosons and 
fermions. By looking at the behavior of xm, fermions 
localize more than bosons. Antisymmetry does help lo- 
calization. This fact that may sound counterintuitive can 
be understood (see the Supplementary Information) by 
looking at the sign of the interference term on the vari- 
ance of Xm- An opposite behavior is observed for the 
distribution of the relative distance R = \j — k\ between 
the two particles. Because of the Pauli exclusion prin- 
ciple, the average distance between the particles in the 
fermionic case is larger than in the bosonic case (where 
they tend to bunch). This is exactly what we observe by 
looking at the variance of the distribution of the relative 
distances (Fig. pp,d), which follows from the probability 



distribution P{R). The different experimental distribu- 
tions of the distance between the particles, obtained in 
both the ordered and static-disordered case, and the dif- 
ferent behaviour of bosons and fermions are shown in Fig. 
St-h. 

So far we discussed the case of static disorder, however 
different types of disorder affect differently the dynamics 
of the entangled pairs. Since our technology is capable of 
implementing arbitrary phase maps in the QWs, a lattice 
with dynamic disorder was produced applying the same 
phase shifts to MZs belonging to a fixed step of the walk 
{4'm,n = </>n7 Vtti), as in Figure 3b of the Supplementary 
Information. Experimental single- and two-photon out- 
put distributions are summarized in Figure [5^,b for the 
bosonic and fermionic case separately. The effect of a 
fully space-correlated dynamic randomness simulates the 
effect of an external classical environment. In this case 
one can show that the limiting distribution is a binomial 
centered in the middle of the spatial axis and with width 
growing linearly with the square root of the number of 
steps. Thus the system undergoes a diffusion process in 
which the propagation becomes equivalent to a purely 
classical random walk. The data of Figure [5^-c indeed 
confirm that in the case of dynamic disorder the walkers 
spread more easily on the BS tree with respect to the 
static case, but spread less with respect to the ballistic 
diffusion of the ordered case. 

The last scenario we considered is the case of both 
space- and time-dependent disorder, which we will call 
fluctuating disorder, realized with random phase shifts 
4^m,n over the entire MZ network (see Fig. 3c of the Sup- 
plementary Information) . In this configuration the diffu- 
sion process leads to a speckle pattern for the two- walkers 
wavefunction (Figure Isli-f). This case shows that the in- 
teraction with an external classical environment quenches 
the localization effect that would be induced by a lattice 
with static disorder. 

Let us note that the experiments are performed on a 
single phase map realization of each disorder. Although 
there are still features that are linked to the particular 
choice of the (randomly picked) phase maps, the number 
of beam splitters is large enough to allow the clear obser- 
vation of the differences between ballistic, diffusive and 
localized regimes. 

We reported on the experimental realization of a quan- 
tum simulator based on discrete quantum walks of entan- 
gled particles in integrated photonic circuits. By prop- 
erly engineering the phase shifts at the output ports of 
the BSs and by changing the number of QW steps, we 
were able to follow in real time the evolution towards An- 
derson localization. The symmetry of the total wavefunc- 
tion (Fermi- or Bose-like) clearly affects the localization 
properties. Fermi statistics helps localization despite the 
fermions tendency to antibunching. The quantum simu- 
lation we performed will help to ascertain the efficiency 
of quantum algorithms with entangled particles on real- 
istic quantum walks. The capability of our technology to 
implement arbitrary phase maps in QWs paves the way 
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to the experimental quantum simulation of the quan- 
tum dynamics of multi-particle correlated systems and 
its ramifications towards the implementation of realistic 
universal quantum computation with quantum walks. 
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FIG. 5: Experimental results: single- and two-particle dis- 
tributions for a 6-step QW in presence of dynamic (a-c) 
and fluctuating (d-f) disorder. Single particle distributions 
have been computed by tracing out the position of one of 
the particles (summing over the columns of the probability 
distribution matrices) The similarities for bosons, fermions 
and single particle with the expected distributions are re- 
spectively SEos = 0-871 ± 0.004, Sp^r = 0.802 ± 0.006 and 
S^i^g = 0.975 =b 0.003 for the QW circuit with dynamic dis- 



order and Sbo 






0.921 ± 0.004, S§^, 



0.852 ± 0.003 and 



0.991 ± 0.002 for the chip with fluctuating disorder. 
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the experimental approach for the simulation of the An- 
derson localization. A.C., R.O., and R.R. fabricated the 
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paper. 
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